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£3 ■ Abstract. Recent developments in nuclear theory allow us to make a con- 

£h . nection between quantum chromo dynamics (QCD) and low-energy nuclear 

physics. First, chiral effective field theory (xEFT) provides a natural hierar- 
chy to define two-nucleon (NN), three- nucleon (NNN), and even four-nucleon 
. interactions. Second, ab initio methods have been developed capable to test 

' these interactions for light nuclei. In this contribution, we discuss ab initio 

no-core shell model (NCSM) calculations for s-shell and p-shell nuclei with 
NN and NNN interactions derived within %EFT. 

<N 

1 Introduction 
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The major outstanding problem in nuclear physics is to calculate properties 



5_i ■ of finite nuclei starting from the basic interactions among nucleons. There are 



two aspects to this problem. First, the basic interactions among nucleons are 
complicated, they are not uniquely defined and there is evidence that more than 
just two-nucleon forces are important. Second, the nuclear many-body problem 
is very difficult to solve. This is a direct consequence of the complex nature of 
the inter-nucleon interactions. 

Interactions among nucleons are governed by QCD. In the low-energy regime 
relevant to nuclear structure, QCD is non-perturbative, and, therefore, hard 
to solve. Thus, theory has been forced to resort to models for the interaction. 
New theoretical developments, however, allow us connect QCD with low-energy 
nuclear physics through promising bridge of xEFT [1] . 

2 Ab initio no-core shell model 

In the ab initio NCSM, we consider a system of A point-like non-relativistic nu- 
cleons that interact by realistic NN or NN+NNN interactions. Unlike in standard 
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shell model calculations, in the NCSM there is no inert core, all the nucleons are 
considered active. Therefore the "no-core" in the name of the approach. Besides 
the employment of realistic NN or NN+NNN interactions, two other major fea- 
tures characterize the NCSM: i) the use of an harmonic oscillator (HO) basis 
truncated by a chosen maximal HO excitation energy N mau ^hf2 relative to the 
unperturbed ground state of the A-nucleon system; ii) the use of effective in- 
teractions. The reason behind the choice of the HO basis is the fact that this is 
the only basis that allows to use single-nucleon coordinates and consequently the 
second-quantization representation without violating the translational invariance 
of the system. The powerful techniques based on the second quantization and 
developed for standard shell model calculations can then be utilized. Therefore 
the "shell model" in the name of the approach. As a downside, one has to face 
the consequences of the incorrect asymptotic behavior of the HO basis. The sec- 
ond feature comes as a consequence of the basis truncation. In order to speed 
up convergence with the basis enlargement, we construct an effective interaction 
from the original realistic NN or NN+NNN potentials by means of a unitary 
transformation. The effective interaction depends on the basis truncation and by 
construction recovers the original realistic NN or NN+NNN interaction as the 
size of the basis approaches infinity. In principle, one can also perform calcu- 
lations with the unmodified, "bare" , original interactions. Such calculations are 
then variational with respect to the HO frequency Q and the basis truncation 
parameter iV max . 

3 Light nuclei from chiral EFT interactions 

Currently the most promising approach to the construction of accurate inter- 
nucleon forces from QCD is the ;\EFT. The %EFT predicts, along with the NN 
interaction at the leading order, an NNN interaction starting at the 3rd order 
(next-to-next-to-leading order or N 2 LO) [HE]) an d even an NNNN interaction 
starting at the 4th order (N 3 LO) [3]. The details of QCD dynamics are contained 
in parameters, low-energy constants (LECs), not fixed by the symmetry. These 
parameters can be constrained by experiment. A crucial feature of ;\EFT is the 
consistency between the NN, NNN and NNNN parts. As a consequence, at N 2 LO 
and N 3 LO, except for two LECs, assigned to two NNN diagrams, the potential 
is fully constrained by the parameters defining the NN interaction. 

We adopt the potentials of the %EFT at the orders presently available, the 
NN at N 3 LO of Ref. [I] and the NNN interaction at N 2 LO [2]. We use ah initio 
NCSM calculations in two ways. One of them is the determination of the LECs 
assigned to two NNN diagrams, cjj and ce [5] that must be determined in A > 3 
systems. Co (ce) is the strength of the NN-n-N (NNN) contact term. The other 
is testing predictions of the chiral NN and NNN interactions for light nuclei. 

It is important to note that our NCSM results through A = 4 are fully 
converged in that they are independent of the -/V max cutoff and the HQ HO 
energy. This is demonstrated in Fig. [TJ where convergence of the 4 He ground-state 
energy using xEFT interactions with and without the NNN terms is shown. For 
heavier systems, we characterize the approach to convergence by the dependence 
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Figure 1. 4 He ground-state energy depen- 
dence on the size of the basis. The HO frequen- 
cies of hO = 28 and 36 MeV were employed. 
Results with (thick lines) and without (thin 
lines) the NNN interaction are shown. The 
full lines correspond to calculations with three- 
body effective interaction, the dashed lines to 
calculations with the bare interaction. 



Figure 2. Relations between cd and ce for 
which the binding energy of 3 H (8.482 MeV) 
and 3 He (7.7f8 MeV) are reproduced, (a) 4 He 
ground-state energy along the averaged curve, 
(b) 4 He charge radius r c along the averaged 
curve. Dotted lines represent the r c uncer- 
tainty due to the uncertainties in the proton 
charge radius. 



of results on iV max and HQ. 

Fig. [2] shows the trajectories of the two LECs cd — ce that are determined 
from fitting the binding energies of the A = 3 systems. Separate curves are 
shown for 3 H and 3 He fits, as well as their average. There are two points where 
the binding of 4 He is reproduced exactly. We observe, however, that in the whole 
investigated range of cd — ce, the calculated 4 He binding energy is within a few 
hundred keV of experiment. Consequently, the determination of the LECs in this 
way is likely not very stringent. We therefore investigate the sensitivity of the 
p-shell nuclear properties to the choice of the cd — ce LECs. First, we maintain 
the A = 3 binding energy constraint. Second, we limit ourselves to the ce> values 
in the vicinity of the point cjj ~ 1 since the values close to the point co ~ 10 
overestimate the 4 He radius. 

While most of p-shell nuclear properties, e.g. excitation spectra, are not very 
sensitive to variations of cjj in the vicinity of the co ~ 1 point, we were able 
to identify several observables that do demonstrate strong dependence on c^. 
For example, the 6 Li quadrupole moment that changes sign depending on the 
choice of cp. In Fig. [3l we display the ratio of the B(E2) transitions from the 10 B 
ground state to the first and the second l + state. This ratio changes by several 
orders of magnitude depending on the co variation. This is due to the fact that 
the structure of the two l + states is exchanged depending on cr> From Figs. [2] 
and El we can see that for cd < —2 the 4 He radius underestimate experiment 
while for cd > the lowest two 1 + states of 10 B are reversed. We therefore select 
cd = — 1 for our further investigation. 

We present in Fig. 2] the excitation spectra of 10 B as a function of -/V max 
for the chiral NN interaction alone. The convergence with increasing iV max is 
quite reasonable for the low-lying states. Similar convergence rates are obtained 
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Figure 3. Dependence of the B 
B(E2;3^0 -» 1^0)/B(E2;3f -> 1^0) 
ratio on cd, with ce constrained by the A = 3 
binding energy fit, for different basis sizes. 
The HO frequency of hfl = 14 MeV was 
employed. 
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Figure 4. B excitation spectra as function 
of the basis space size iV max at hQ — 15 MeV 
using the chiral NN interaction compared to 
experiment. The isospin of the states not ex- 
plicitly depicted is T=0. 



for other p— shell nuclei that we investigate. A remarkable feature of the 10 B 
results is the observation that the chiral NN interaction alone predicts incorrect 
ground-state spin of 10 B. The experimental value is 3 + 0, while the calculated one 
is l + 0. On the other hand, once we also include the chiral NNN interaction in the 
Hamiltonian, which is actually required by the ^EFT, the correct ground-state 
spin is predicted. Further, once we select the cd value as discussed above, i.e. 
cd = — 1, we also obtain the two lowest l + states in the experimental order. 

We display in Fig. [5] the natural parity excitation spectra of four nuclei in 
the middle of the p— shell with both the NN and the NN+NNN interactions 
from xEFT. The results shown are obtained in the largest basis spaces achieved 
to date for these nuclei with the NNN interactions, N max = 6 (6hf2). Overall, 
the NNN interaction contributes significantly to improve theory in comparison 
with experiment. This is especially well-demonstrated in the odd mass nuclei for 
the lowest few excited states. The case of the ground state spin of 10 B and its 
sensitivity to the presence of the NNN interaction discussed also in Fig. HI is 
clearly evident. We note that the 10 B results shown in Fig. U] were obtained with 
the HO frequency of HQ = 15 MeV, while those in Fig. [5] with hfi = 14 MeV. A 
weak HO frequency dependence of the -/V max = 6 results is apparent. 

Further details on calculations presented in this section were published in 
Ref. [6]. These calculations demonstrate that the chiral NNN interaction makes 
substantial contributions to improving the spectra and other observables. How- 
ever, there is room for further improvement in comparison with experiment. We 
used a strength of the 27r-exchange piece of the NNN interaction, which is consis- 
tent with the NN interaction that we employed (i.e. from Ref. [3]). This strength 
is somewhat uncertain (see e.g. Ref. [5]). Therefore, it will be important to study 
the sensitivity of our results with respect to this strength. Further on, it will be 
interesting to incorporate sub-leading NNN interactions and also four-nucleon 
interactions, which are also order N 3 LO [3]. Finally, it will be useful to extend 
the basis spaces to N max = 8 (8hQ) for A > 6 to further improve convergence. 
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Figure 5. States dominated by p-shell configurations for 10 B, 11 B, 12 C, and 13 C calculated at 
A^max = 6 using hQ = 15 MeV (14 MeV for 10 B). Most of the eigenstates are isospin T=0 or 
1/2, the isospin label is explicitly shown only for states with T=l or 3/2. The excitation energy 
scales are in MeV. 

4 Beyond nuclear structure with chiral EFT interactions 

A shortcoming of the ab initio NCSM is its incorrect description of long-range 
correlations and its lack of coupling to continuum due to the expansion of the 
eigenstates in a finite HO basis. If we want to build upon the ab initio NCSM 
to microscopically describe nuclear reactions we can proceed in two ways. First, 
we can rely on techniques such as the Lorentz Integral Transform that reduce 
the continuum problem to a bound-state-like problem. Using this approach, we 
performed 4 He photo-disintegration cross section calculations using the %EFT 
interactions [7] . We demonstrated a sizeable effect of the chiral NNN interaction 
on the cross section. Second, we can augment the ab initio NCSM basis by ex- 
plicitly including cluster states and solve for their relative motion while imposing 
the proper boundary conditions. This approach, applicable to a wide range of re- 
actions as well as to weakly bound states, is in the spirit of the resonating group 
method (RGM) [8]. In our approach, we use the ab initio NCSM wave functions 
for the clusters involved and effective interactions derived from realistic NN (and 
eventually also from NNN) potentials. As an example, a converged calculation 
of the n+ 4 He 2 Si/2 phase shift using the x^FT NN interaction is presented in 
Fig. [6j More details on this approach are given in Ref. [DJ. 

5 Conclusions 

The ab initio NCSM evolved into a powerful many-body technique. Presently, 
it is the only method capable to use interactions derived within the x^FT for 
systems of more than four nucleons. Among its successes is the demonstration 
of the importance of the NNN interaction for nuclear structure. Applications 
to nuclear reactions with a proper treatment of long-range properties are under 
development. Extension to heavier nuclei is achieved through the importance- 
truncated NCSM [ID] . Within this approach, ab initio calculations for nuclei as 




heavy as Ca become possible. 

Acknowledgement. This work performed under the auspices of the U.S. Department of Energy 
by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. Support 
from U.S. DOE/SC/NP (Work Proposal Number SCW0498) and the Department of Energy 
under Grant DE-FC02-07ER41457 is acknowledged. 



References 

1. S. Weinberg, Physica 96A, 327 (1979); Phys. Lett. B 251, 288 (1990); Nucl. 
Phys. B363, 3 (1991); J. Gasser et al, Ann. of Phys. 158, 142 (1984); Nucl. 
Phys. B250, 465 (1985). 

2. U. van Kolck, Phys. Rev. C 49, 2932 (1994); E. Epelbaum, A. Nogga, W. 
Glockle, H. Kamada,Ulf-G. Meissner and H. Witala , Phys. Rev. C 66, 064001 
(2002). 

3. E. Epelbaum, Phys. Lett. B 639, 456 (2006). 

4. D. R. Entem and R. Machleidt , Phys. Rev. C 68, 041001(R) (2003). 

5. A. Nogga, P. Navratil, B. R. Barrett and J. P. Vary Phys. Rev. C 73, 064002 
(2006). 

6. P. Navratil and V. G. Gueorguiev and J. P. Vary, W. E. Ormand and A. 
Nogga, Phys. Rev. Lett. 99, 042501 (2007); [m!cLth/0701038[ 

7. S. Quaglioni and P. Navratil, Phys. Lett. B 652, 370 (2007). 

8. Y. C. Tang, M. LeMere and D. R. Thompson, Phys. Rep. 47 (1978) 167; K. 
Langanke and H. Friedrich, Advances in Nuclear Physics, chapter 4., Plenum, 
New York, 1987. 

9. S. Quaglioni and P. Navratil, larXiv:0712.0855l 

10. R. Roth and P. Navratil, Phys. Rev. Lett 99, 092501 (2007). 



